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Gain in current-driven semiconductor heterostructure devices is calculated within the theory of 
nonequilibrium Green functions. In order to treat the nonequilibrium distribution self-consistently 
the full two-time structure of the theory is employed without relying on any sort of Kadanoff-Baym 
Ansatz. The results are independent of the choice of the electromagnetic field if the variation of the 
self-energy is taken into account. Excellent quantitative agreement is obtained with the experimental 
gain spectrum of a quantum cascade laser. Calculations for semiconductor superlattices show that 
the simple 2-time miniband transport model gives reliable results for large miniband widths at room 
temperature. 
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I. INTRODUCTION 



The prospect of a semiconductor laser in the infrared 
and THz region has been one of the key reasons for the 
development and study of semiconductor heterostructure 
elements, since the first proposal of semiconductor super- 
lattices in 1970 13]. A possible gain mechanism may be 
based on two different ideas: (i) At certain electrical fields 
resonant tunneling between different subbands can lead 
to population inversion associated with gain at the tran- 
sition energy [2]. This idea was realized in the quantum 
cascade laser [3], which has become an important device 
in the infrared region. Lasing in the THz region has 
been demonstrated very recently as well [Q . For further 
details, see the review 

§■ (ii) The occurrence of nega- 
tive differential conductivity in superlattices gives raise 
to gain in the low frequency range extending up to fre- 
quencies of the order of the Bloch frequency Q . Despite 
a strong effort by several groups all over the world, this 
idea is still not realized. The difficulty in its realization is 
attributed to the instability of the operating state lead- 
ing to domain formation jj, which may be circumvented 
by different strategies [§, f|, 0. For further references 
and a detailed discussion sec Ref. [TJ. 



These concepts for gain are commonly described in dif- 
ferent ways: (i) An external electromagnetic field causes 
transitions between different levels, where the transition 
rate R is evaluated by Fermi's golden rule. If stimulated 
emission dominates stimulated absorption (i.e., for popu- 
lation inversion), gain occurs at the transition frequency. 
For a pair of states a, (3 with energies E ai Ep, occupation 
probabilities f a , fp, and a dipolc matrix element z a p, the 
contribution to the material gain (i.e., increase of light 



intensity per length) is given by 

1 K a 



G m (w) = 



■stim. em. pstim. abs. 

V Photon flux per area 

7TW \eZa/3 



(1) 



e r ceoV 



8(E a -Ep-hw){f a - f p ). 



where V is the normalization volume, e < is the elec- 
tron charge, c is the vacuum speed of light, and e,. the 
background dielectric constant. 

(ii) Transport theory in the presence of alternating 
electric fields provides the complex dynamical conduc- 
tance a(ui). Standard electrodynamics (see, e.g., section 
7.5 of Ref. Jl^l ) gives the material gain (the negative 
absorption coefficient a, which is assumed to be small 
compared to the wave vector here) 



3?{a( W )} 



(2) 
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Thus, gain is equivalent to a negative differential conduc- 
tivity in the respective frequency range. This concept has 
been frequently applied to supcrlattice transport. (For 
superlattices Eq. ([!]) gives zero gain in the basis of Wan- 
nier Stark states, as the translational symmetry gives 
identical f a for all Wannier-Stark levels.) Here the cur- 
rent is either evaluated from a semiclassical miniband 
transport model ||, ^, ^, |lC|_ or by sequential tunnel- 
ing between different layers i|| [D| [uj. While these ap- 
proaches use a macroscopic current density in Eq. (Q), the 
consideration of polarization currents allows for a deriva- 
tion of Eq. (|I]), see, e.g., Ref. |l5[ 

For complicated semiconductor heterostructures, with 
a variety of different tunneling and optical transitions, 
gain may result both from optical transitions in the spirit 
of (i) and macroscopic currents in the spirit of (ii). The 
aim of this paper is to show the feasibility of a general 
approach based on a quantum transport theory which 
treats both mechanisms on a equal footing. 

The method of nonequilibrium Green functions (T(| 
is used here, which has been frequently applied to light 
emission in semiconductors (see, e.g., Refs. O, nM H9, EG 



2 



and references given therein) . While these works consider 
transitions between the conduction and valence band, the 
focus in this paper is on transitions within the conduction 
band in current-driven heterostructure devices. The cal- 
culations are performed in two steps: First the transport 
problem is solved by evaluating a self-consistent station- 
ary solution of the quantum kinetic equations for a given 
applied bias. This provides us with the Green functions 
describing the electronic state far from equilibrium. In a 
second step, an additional weak radiation field is taken 
into account. The time-dependent kinetic equations are 
linearized around the stationary nonequilibrium state in 
oder to study the linear response of the current-driven 
system. The formulation used here employs the full two- 
time structure of the theory without any sort of (gen- 
eralized) Kadanoff-Baym ansatz The capability of 
the approach is demonstrated by calculations of the gain 
spectra in a quantum cascade laser and a superlattice. 



II. THEORY 

We describe a general system by a set of orthonormal 
states labeled a with energies E a and write the Hamil- 
tonian in the form: 



a a/3 



(3) 



where a a and a^ a are electron annihilation and creation 
operators in the state a. Uapit) represents the matrix el- 
ements of the kinetic energy and the potential part of the 
Hamiltonian H. In the following, U(t) will be split into 
a constant part U and a small time-dependent pertur- 
bation 5U (t) describing the interaction with a radiation 
field. Finally -ff sca tt contains the part of the Hamiltonian 
which will be solved perturbatively (such as impurity, 
phonon, or electron-electron scattering matrix elements). 

Within the theory of nonequilibrium Green functions 
]16| , |l8| the key quantities are the correlation function (or 
'lesser' Green function) 

G< >aa (t 1 ,«2) = i(aL(*a) a « l (*i)) W 
and the retarded/advanced Green functions 

G^(t 1 ,t 2 ) = T ie[±(ti-t2)] 

I . , \ (5) 

x {a ai (ii)aT 2 (t 2 ) + al 2 (t 2 )a ai {tx)J, 

respectively, where the Heisenberg picture is used. 

First consider a stationary state with a time indepen- 
dent matrix U{t) — U, neglecting the radiation field. In 
this case all functions depend only on the time differ- 
ence ti — t 2 [these stationary state functions are labeled 
by a tilde, e.g., G{t\ — t 2 ), in the following], and it is 
convenient to work in Fourier space defined by 



G aua2 (E) = ± J dte im / h G ai , a2 (t) 



(6) 



Then Eqs. (A3 A4) of the appendix become the matrix 
equations 



(E - E ai ) G™* Q2 (E) - ^ U ai ,p Gff a2 (E) 



(J) 



(8) 



These have to be solved self-consistently together with 
the equations for the self-energies which are functionals 
of the Green functions 



t(E) =J r E {G Iet (E'),G < (E')\ 



(9) 



This standard approach allows for a self-consistent eval- 
uation of the Green functions G in a nonequilibrium sit- 
uation caused by an applied bias (contained in U). De- 
tails, such as the specific form of the functionals in the 
self-consistent Born approximation, can be found in, e.g., 

The current density (in the z-direction) can be evalu- 
ated from the expectation value of the momentum oper- 
ator divided by the mass: 



e i 



Vh 



{[H,z]) 



' ' " ( |0] 

= — 22 ^af}G^ a (t, t) + Scattering currents 



a/3 



where W a p = ^Z^U^z^p — z a ~ ( U 7 p). The scattering cur- 
rents result from the scattering part of the Hamiltonian 
and can be expressed in terms of self-energies. Details 
will be given elsewhere. 

Now we consider the influence of an additional time- 
dependent potential 



5U(t) = J ^SV{uj)e- 



(11) 



in the Hamiltonian. (Bold capital symbols denote matri- 
ces in the state indices a,/3). We treat the change of the 
system in linear response, and set 



G(t 1 ,t2) = G(t l -t 2 ) + SG(t 1 ,t 2 ). 



(12) 



The chance in the self-energies is described within the 
linearization of the functional (here in the time domain) 

£(ti, t 2 ) =Ft{G + 5G} « E(ti - ta) + <JE(ii, t 2 ) . 

(13) 

Such a decomposition has been used in Ref. [l?] for 
Hartree-Fock self-energies, e.g.. In this case, the time- 
dependence of the self-energies allows for a reduction to 
density matrix equations by setting t\ =t 2 . In contrast, 
in our case, where scattering effects are considered, 5G 
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exhibits an explicit time dependence in both arguments. 
This two time structure is fully taken into account here. 
Therefore we apply the Fourier decomposition in both 
times via 

8G( tl ,t 2 )= [ [ ^SG(u;,E)e- iE ^y h 

J 2ir J z7r 

(14) 



The same decomposition is used for <5£. It is shown in 
appendix ^ that 5G is determined within linear response 
by the following equations: 



5G ret/adv (u;,£0 



G iet/adv (E + hu) SU(u) + <5£ rot/adv (w, E) G ret/adv (.E) 

G rct (E + haj)S\J{uj)G < (E) + G<(E + hw)5U(uj)G adv (E) 

G TCt {E + ^)JS rct (w, E)G<(E) + G rct (E + ftw)ffi<(w, ^)G adv (^) 

G<(E + hw)Si: adv (cj, E)G adv (E) 



(15) 



(16) 



The changes in self-energy <5£ are functionals of 5G, 
which have to be evaluated self-consistently with SG. 
The derivation of the functionals from the linearization 
of Eq. ( |l3| ) is straightforward and can be performed for 
arbitrary self-energies. The situation is particularly sim- 
ple for the self-consistent Born approximation, where the 
functional T is linear in G, and one finds: STi(uj,E) = 
F E {SG Ict (uj,E'),SG < (uj,E')} with the same functional 
as Eq. @. (Note, that c5G adv (w, E') = [SG Ict (uj, E')] f 
does not hold for the time-dependent quantities. Thus, 
the advanced quantities should be calculated explicitely) 

The terms related to SH correspond to the ladder cor- 
rections in the evaluation of diagrams for the Kubo for- 
mula in equilibrium. It is a general advantage of lin- 
ear response within nonequilibrium Green functions that 
these terms are obtained directly, see also |l8|, p^ |. 

Now we consider the linear response to an external 
radiation field where the electric field points in the z di- 
rection. This gives additional terms in the Hamiltonian 
H = (p — eA) 2 /2m + V(r) + eip(r), where A,<p are the 
electromagnetic potentials. Neglecting quadratic terms 
in the radiation field and terms containing k (i.e., assum- 
ing that the wavelength is large compared to the size of 
the active region), the following perturbation potentials 
are obtained (see Appendix : In the Coulomb gauge 



6U a , p {u) = ^F{u)(^- 

1U1 \ TO 



a, (3 



eF{uj) 



a,0 



(17) 



neglecting the scattering part of H; in the Lorentz gauge 



5U a ,p(u) = -eF(cj)z a ,fi ■ 



(18) 



change in Eq. ( pL 0| ) 
W.M = ^J Tr{ [SUM Z Z 6V(cj)\ 

■ G<(E) +W • 5G < (w,E)\. 

(19) 



Finally, the change in current density is given by the 



The scattering current in Eq. (10) has not been included 
here. With these ingredients, one obtains the complex 
conductivity o-(u>) = 5 J (u>) / F '(cj) , as well as the gain 
coefficient via Eq. (0). 



III. RESULTS 

In the following we consider two different structures: 
(i) The quantum-cascade-laser structure of Ref. |||. (ii) 
The superlattice used in Ref. ^ojwith 2 nm Alo.3Gao.7As 
barriers and 8 nm GaAs wells and a doping of n = 
10 16 /cm 3 . In both cases, we use a basis set consisting of 
products of Wannicr functions (in the growth direction 
z) and plane waves with wave vector k (in the (x, y)- 
plane) , and obtain U a p from nominal sample parameters 
[jllj . The self-consistent solution for the electrostatic po- 
tential is included as well. 

For the scattering Hamiltonian we include interface- 
roughness scattering and phonon scattering (including 
both optical phonons and a second phonon with low en- 
ergy to mimic acoustic phonons). The self-energies are 
evaluated in the self-consistent Born approximation. Fur- 
thermore, we use momentum-independent scattering ma- 
trix elements (evaluated for a typical momentum trans- 
fer) and diagonal self-energies. In this approximation the 
self-energies are not k-dependent, which significantly re- 
duces the numerical effort. Results and further details 
for superlattices |0, and quantum cascade lasers p7j 
have been given previously. The current-field characteris- 
tics are displayed in Fig. |l|(a,b) for both structures, which 
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FIG. 1: Current density versus potential drop per period for 
the quantum cascade laser of Ref. ^ (a) and the superlat- 
tice structure of Ref. (b). The full line gives the result 
of the quantum transport model (QT), while the dashed line 
is evaluated with the simple 2-time miniband model (MB), 
Eq. (Ek. 



agree reasonably well with the respective experimental 
data. Note that domain formation in the region of neg- 
ative differential conductivity modifies the behavior of 
the experimental curve for the superlattice. Such effects 
are beyond the scope of the general approach discussed 
here, but can be easily studied using simpler transport 
models jy], |2^]. Similar results for the current-field 
relation of this quantum cascade laser have been reported 
in Ref. |o[ 

For the superlattice structure, the simple 2-time mini- 
band transport model |3lj gives 



Jmb(F) = en 



dAI 1 (A/2k B T) eFdh/T e 



2h I (A/2k B T) (eFd) 2 + ft 2 /r e r„ 



(20) 



for a non-degenerate electron gas. Here F is the electri- 
cal field, A the miniband width and d the period of the 
superlattice. r e = 0.17 ps and r m = 0.113 ps represent 
fitted, phenomenological energy and momentum scatter- 
ing times. Ij(x) are the modified Bessel functions. In 
Fig. 0(b) this model shows reasonable agreement with the 
full quantum transport calculation for |eFd| < A/2 = 11 
meV, the field range where miniband transport holds true 
p2| . (For lower temperatures the shape of the current- 
field relation becomes more complicated both in the semi- 
classical Boltzmann and the quantum transport model 
p6| so that the simple 2-time model cannot fit the data.) 

In the following the response to an external radiation 
field is studied: First let us use the Coulomb gauge and 
neglect the terms with ST, in Eq. (|l6|). The gain spec- 
trum for the quantum cascade laser of Ref. 25 is displayed 
in Fig. 0. Gain sets in for current densities above 0.8 
kA/cm and increases with current. The peak gain co- 
efficient is 57 cm" 1 at 6.5 kA/cm 2 and a photon energy 
of 130 meV, in excellent agreement with the findings of 
Ref. 25 (estimated losses divided by confinement factor 
63 cm" 1 and threshold current J t h=7.2 kA/cm 2 , lasing 
at 131 meV). The width of the gain spectrum agrees well 
with the findings of Ref. |3^. 

Results for the superlattice are given in Fig. ||. The 
real part of the conductivity essentially follows the result 
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FIG. 2: Gain spectrum for the quantum cascade laser of 
Ref. ^5] evaluated from the quantum transport model with 
Coulomb gauge and neglecting self-energy corrections. 
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FIG. 3: Real part of the dynamical conductance for the super- 
lattice structure of Ref. hoi. Full and dashed line: Result from 
the quantum transport model with Coulomb gauge and ne- 
glecting self-energy corrections for Fd = 5 mV and 20 mV, 
respectively. Dotted and dash-dotted line: Corresponding re- 
sults from the 2-time miniband model [HI. 



from the simple 2-time miniband model 
Jmb(F) 1 - uj B T rn T e - k 



a(u>) 



F 



(-1 



i0 2 )T m T e + 1 - ko(T m + T e ) 



(21) 



with uj b — eFd/h. Thus, this simple model gives, at 
least for wide minibands and high temperatures, reliable 
results. 

Now we want to study the relevance of the terms 5T, 
in Eq. (|l^), which have been neglected so far. Fig. |] 
shows that their inclusion barely changes the result in the 
Coulomb gauge (the full and dashed lines are almost in- 
distinguishable). Therefore these self-energy corrections 
may be neglected (at least for diagonal self-energies as- 
sumed here) allowing for a significant reduction of the 
numerical effort. In contrast, in the case of the Lorentz 
gauge |Q, the terms evaluated with self-energy correc- 
tions (dash-dotted line) differ substantially from the sim- 
pler version without these terms (dotted line). The differ- 
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FIG. 4: Comparison of different choices of the gauge and self- 
energy corrections ST, for (a) the superlattice at Fd — 20 and 
(b) the quantum cascade laser at Fd = 260 mV. Full line: 
Coulomb gauge neglecting <5£; Dotted line: Lorentz gauge 
neglecting 5E; Dashed line: Coulomb gauge with 8T, (falls 
together with the full line); Dash-dotted line: Lorentz gauge 
with <5E. 



of the choice of the gauge. 



IV. CONCLUSION 

A formulation of the theory of nonequilibrium Green 
functions has been given which allows for a quantitative 
description of current-induced gain in semiconductor het- 
erostructures. Quantitative agreement with the experi- 
ment for the gain spectrum of a quantum cascade laser 
was obtained. The standard 2-time miniband model for 
the superlattice could be verified for a large miniband 
width and under room temperature operation. The nu- 
merical results are not sensitive to the choice of the gauge 
if the variations of self-energy terms are taken into ac- 
count. These terms are of particular importance for the 
Lorentz gauge, while their contribution is small in the 
Coulomb gauge in the calculations presented here. 



ence is particularly strong in the case of the superlattice, 
where the diagonal terms in the perturbation potential 
( |l8| ) are essential for the transition. Furthermore, the 
results of Coulomb gauge and Lorentz gauge are almost 
identical if the self-energy corrections are included. This 
demonstrates that the inclusion of these terms is essential 
to obtain a consistent theory, which must be independent 
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APPENDIX A: DERIVATION OF THE EQUATIONS FOR THE STATIONARY STATE 

The Green functions are determined by the Dyson Equation (see, e.g., chapter 5 of Ref. |l^) 

' P J 

= h8(tx - t 2 )S aua2 (Al) 

and the Keldysh relation, see also Appendix [b|. 

0,0' 



(A2) 



Now we consider the stationary state, where G(t\,t 2 ) — G{t\ — t 2 )- Eqs. (Al A2) become: 

n^gr (*) - E v^&g!?® - E / ^-KlTt - *)&?J?W) - 



g<, Q2 w-E/x/x g ^ 

0,0' 



(t-t')X<(t'-t")Gpv(t") 



(A3) 
(A4) 



and their Fourier-transformation (||) provides us with Eqs. (0J|). 



APPENDIX B: A REMARK ON EQ. (5.11) OF REF. ^| 



Eq. (A2) is a particular solution of the inhomogeneous differential equations (5.3) of Ref. [l8|, which are given more 
explicitly in Eqs. (106-107) of Ref. O in the notation used here. The Keldysh relation (5.11) derived in Ref. Il8| 
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contains an additional term: 

G< dd (t 1 ,t 2 ) = J^J^ hS(t 1 -t b )l + G Ict (t u t b )U(t b ) + J ^G ret (ti, t )£ ret (t 0) t h ) 

xG<(t b ,t c ) h5(t c -t 2 )l + U(t c )G lldv (t c ,t 2 ) + J ^^ ct (t c ,t d )G* dv (t d ,t 2 ) 



(Bl) 



This term constitutes a homogeneous solution of the inhomogeneous differential equations (for fixed E and G? ret / adv ) 
which guarantees that G < — ► Gq for t\ 1 t 2 - * — oo. This satisfies a general assumption underlying the perturbation 
expansion of the S'-matrix. In the following we show that G^ dd vanishes for finite times t\,t 2 . 
Using Eqs. (108,109) of Ref. O one obtains 



G< dd (h,h) -J^J ^G^(t llta ) (-i^ - E ) G<(i Q ,i c ) (ih 



d_ 

dt P . 



-E)G adv (t c ,t 2 ) 



7x/ t Grct ^'^«- E ) G o^^«- E ) Gadv 



(tc,t 2 ) 



(B2) 



dt c 







G Iet (ti,t a )G$(t a ,t c ) ( ih— - E ) G a(iv (t c ,t 2 ) 



where d means that the derivative operates to the left. In the second line, partial integration has been used. As 

(ih--^ E^j Gq (t a ,t c ) = 0, the first term vanishes. Furthermore, at the upper bound G ret (ii,oo) = holds for finite 

times t\. Applying the same manipulations for t c , we find 



G< d d(^2) =i J ^G ret (i 1 ,-oo)G Q < (-oo,t c ) ^A-E^ G adv (i c ,t 2 ) 
=G rct (i 1 , -oo)G<(-oo, -oo)G adv (-oo, t 2 ) . 



(B3) 



For a real system the retarded and advanced Green functions decay for large time differences and we find G^ dd (t% , t 2 ) — ► 
for finite t\, t 2 . Therefore this additional term is spurious. 

A second, more physical argument relates to the idea, that the stationary state of a real physical system should 
not depend on the initial conditions used in a formal theory. Except for specific situations exhibiting a hysteresis, 
scattering processes will install a stationary state which is independent of initial conditions. Therefore all additional 
homogeneous solutions should reflect transient effects, which vanish if the initial state is chosen at t\ = —00 and 
t 2 = —00. 



APPENDIX C: DERIVATION OF EQS. ([ij,[l6|) 



Setting U(t) = U + SU{t), G(t 1 ,t 2 ) = G(h - t 2 ) + 5G{t u t 2 ) and E(h,t 2 ) = £(ii - t 2 ) + 5T,(h,t 2 ) in Eqs. (|Al].|A|) 
gives in linear order 5U: 



^SGT ua2 (h,t 2 ) J2 U au ^% 2 {t u t 2 ) - W ^E£,0(ti - t)8G)^{tM) = 

£* c W*0<3£ t a,(*i-*2)+E / T^A^GfaQ-**) (CI) 

P J 

This constitutes an inhomogeneous differential equation, where the homogeneous part (left-hand side) is identical 
with the defining equation (A3) for G rot . Thus G ret is the corresponding Green function and the solution is given by: 

5G iet (h,t 2 ) = fjf ^-G rot (ti-*)[<5U(t)fi*(t-* / ) + <yS ret (*,*')]G ret (t'-*2) (C2) 



and Fourier transformation using Eq. ( |l4| ) gives Eq. ([[5]) after some algebra. SG adv is evaluated in the same way with 
ret replaced by adv. 
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The linearization of the Keldysh relation (|A2| ) provides us with 

5G<(t 1 ,t 2 ) = J^-J ^^G rct (t 1 ,t)S < (t-OG adv (t'-i2) + G rc Hii-<)^<(t,tOG adv (t'-f 2 ) 

+ G^ih - t)t<{t - t^SG^it' ,t 2 )] (C3) 

Fourier transformation gives Eq. ( |l6|) after some algebra. 

APPENDIX D: CHOICE OF GAUGES 

We study the linear response to an external radiation held propagating in y-direction 

p(f,t)= [ ^-F(uj)e z e ik ^y- iut B(r,t)= [ ^^F(u)eJ k W y - iut (m) 

J 2lT J Z7T UJ 

Thus we have two reasonable choices of the electromagnetic potentials [Ell : 

Coulomb gauge A(f, t) = f ^¥e/ W9 " w (f>(f,t) = (D2) 

J 2ir iui 

Lorentz gauge A{r,t) = - [ ^^^^ y^^^'"' 4>(r, t) = -F{uj)ze lk( - u})v ~ ibJt (D3) 

J Z7T LU 

In the spirit of a dipole approximation (i.e., assuming that the wavelength is large compared to the size of the active 
region) the terms containing k are neglected, providing the expressions in the perturbation Hamiltonian, where 

quadratic terms cx F(uj) 2 are neglected. 
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